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Abstract 

A generic lattice model for systems containing particles interacting with short-range attraction 
long-range repulsion (SALR) potential that can be solved exactly in one dimension is introduced. 
We assume attraction J\ between the first neighbors and repulsion Ji between the third neighbors. 
The ground state of the model shows existence of two homogeneous phases (gas and liquid) for 
J2/J1 < 1/3- In addition to the homogeneous phases, the third phase with periodically distributed 
clusters appears for J2/J1 > 1/3. Phase diagrams obtained in the self-consistent mean- field ap- 
proximation for a range of values of J2/J1 show very rich behavior, including reentrant melting, and 
coexistence of two periodic phases (one with strong and the other one with weak order) terminated 
at a critical point. We present exact solutions for the equation of state as well as for the correlation 
function for characteristic values of J2I J\- Based on the exact results, for J2/J1 > 1/3 we predict 
pseudo-phase transitions to the ordered cluster phase indicated by a rapid change of density for 
a very narrow range of pressure, and by a very large correlation length for thermodynamic states 
where the periodic phase is stable in mean field. For 1/9 < J2/J1 < 1/3 the correlation function 
decays monotonically below certain temperature, whereas above this temperature exponentially 
damped oscillatory behavior is obtained. Thus, even though macroscopic phase separation is ener- 
getically favored and appears for weak repulsion at T = 0, local spatial inhomogeneities appear for 
finite T. Monte Carlo simulations in canonical ensemble show that specific heat has a maximum 
for low density p that we associate with formation of living clusters, and if the repulsion is strong, 
another maximum for p = 1/2. 
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I. INTRODUCTION 



Most of biologically relevant macromolecules, or particles in soft-matter systems are 
charged and repel each other with screened electrostatic forces jl-4]. On the other hand, 
complex solvents in biological or soft-matter systems may induce effective attraction between 
the macromolecules. Important examples of the effective attraction include depletion forces 
resulting frointhe presence of small objects like nonadsorbing polymers {2, 5, 6], solvophobic 



12J resulting from critical adsorption 



attraction {7, 8] or thermodynamic Casimir forces 
near a critical demixing point in a binary (or multicomponent) solvent. The sum of all inter- 
actions has often the form of short-range attraction long-range repulsion (SALR) potential 
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Despite of the importance of the SALR potential for biological and soft-matter systems 
only preliminary and incomplete results for phase diagrams have been obtained so far 



261 ] . It is clear that the phase diagrams depend on the strength of the repulsion, and for 



strong enough repulsion can be completely different than in simple fluids. In the latter 
systems the phase behavior is determined by the competition between the entropy favoring 
the disordered phase and the energy favoring formation of a spherical droplet of the dense 
phase. Volume fraction of the particles determines only the size of this droplet. In the case 
of the SALR potential there is additional competition between the energy and the chemical 
potential (or the volume fraction) of particles. The SALR potential favors formation of 
spherical clusters whose size is determined by the range of the attraction, and the distance 
between the clusters is determined by the range of the repulsion. The structure minimizing 
the energy is possible only for sufficiently small volume fraction, however. The increase 
of the chemical potential leads to transitions from spherical to cylindrical clusters, next to 
a network, then to layers of par ticles, and finally to inverse structures (occupied regions 
replaced by voids) l2ll . |23| . |24J . |26|. Regions of stability of the disordered fluid and 



different ordered phases, as well as location and order of phase-transitions between them for 
various shapes of the SALR potential are only partially known. 

The difficulty in determination of phase diagrams in experimental studies results from the 
presence of many metastable phases and large characteristic time scales, by which instan- 
taneous rather than average states have been observed 

HQ . The question of the average 
structure in thermodynamic equilibrium may be clarified in future confocal microscopy dy- 
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namical studies. In simulations it is much more difficult to study inhomogeneous structures 
with some kind of periodic ordering and moreover the results are restricted to the chosen 
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271 ] . In theoretical studies the results relay on approxi- 
25]. Validity of MF for the SALR type 



form of interactions 

mations, usually of the mean-field (M F) ty 
potentials was only partially verified 20j, ]23J, and it is not clear which predictions of the 
MF approximation and to what extent are valid on a qualitative or on a semiquantitative 
level. Preliminary results within field-theoretic a ppr oach show that fluctuations may lead to 
substantial modification of the equation of state 



and to a fluctuation- induced first order 



phase transition 
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29], as in magnetic systems with competing interactions [30l . |31 |- We 



should also mention that the accurate liquid theories such as SCOZA 32J are unfortunatelly 
limited to uniform phases, and only a lack of solutions is an indicator of possible phases 
with periodic ordering on the mesoscopic length scale [lj],[l8[ 19]. Further studies including 
comparison between MF type results and simulations, as in Ref. 23] are necessary. It would 
be also desirable to introduce a simple model that could be solved exactly. Exact results 
could serve for testing general theoretical predictions and for verification of the range of 
validity of approximate theories. No simple model that could play the same role as the two 
dimensional Ising model played for simple systems was proposed so far. 

An exactly solvable two-dimensional model for the SALR potential cannot be proposed at 
the moment. To fill this gap at least partially, we introduce in this work a one-dimensional 
(Id) lattice model with attraction J\ between the first and repulsion J2 between the third 
neighbors. Exact solutions can be obtained for the equation of state (EOS) as well as 
for the density and the correlation function for the whole range of Jij J\. There are no 
phase transitions in Id models with finite range of interactions, but pseudo-phase transitions 
may exist for low temperatures. Moreover, it is interesting to compare the exact and MF 
results. It is also of interest to compare thermodynamics and structure of the disordered 
phase in simple fluid and in the SALR-potential systems. The results can give interesting 
information about pretransitional ordering in the latter systems. Effects of the precursors of 
the ordered structures on structure, compressibility and specific heat were already studied 



in approximate theories and in simulations 
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251 ] . Based on the exact results for the 



correlation function and for the EOS, we can verify if the repulsion can lead to significantly 
increased compressibility and to oscillatory decay of correlations on a mesoscopic length 



scale, as found in Ref. 
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23]. 
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In MF phase transitions are present even in one dimension, therefore MF in Id systems 
can shed light on some general features of phase diagrams in two or three dimensions. 

Formation of living (or dynamic) clusters may be associated with large energy fluctua- 
tions, especially when the clusters easily form or dissociate. The specific heat in the SALR 
system was studied in Ref.[l7|]. For fixed densities the maximum of the specific heat at 
some temperature was interpreted as a signature of the transition between the homogeneous 
and inhomogeneous fluids. To calculate the specific heat in our model we perform Monte 
Carlo (MC) simulations in the canonical ensemble. We calculate the specific heat for several 
temperatures as function of density for strong and for weak repulsion. 

The one dimensional model considered here, apart from yielding general information on 
periodic ordering in the SALR-potential systems, describes several physical systems that are 
interesting by their own. In the first place it can represent charged particles in a presence 
of depletant at a three-phase coexistence line, or adsorbed at nanotubes or microtubules. 
Another example is a linear backbone polymer with monomers containing sites binding 
particles or ions that attract or repel each other when bound to first or third neighbors on 
the backbone respectively. Our model can answer the question of spontaneous formation of 
ordered periodic patterns on linear substrates. 

In the next section the model is introduced and the ground state is determined and 
discussed. Sec. 3 is devoted to the MF approximation. Exact results and MC simulations 
are described in sees. 4 and 5 respectively. Sec. 6 contains a short summary and discussion. 

II. THE MODEL AND ITS GROUND STATE 
A. The model 

We consider an open system in equilibrium with a reservoir with temperature T and chem- 
ical potential \x v . The interaction h between the particles and the nanotubes, microtubules 
or binding sites plays analogous role as the chemical potential, and we introduce p = p p + h. 
The particles can occupy lattice sites labeled by x taking integer values, 1 < x < L and 
we assume periodic boundary conditions (PBC), i.e. L + 1 = 1,0 = L. Each microstate is 
described by {p(x)} = (p(l), ...,p(L)), where the occupancy operator p(x) = 1 or p(x) = 
when the site x is occupied or empty respectively. The probability of the microstate {p(x)} 
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is given by 

-/3H[{p(x)}] 

pmm = — = — , (i) 

where S is the normalization constant, and /3 = (fceT) -1 with k B denoting the Boltzmann 
constant. We have introduced for convenience the thermodynamic Hamiltonian containing 
the energy and the chemical potential term, 

H [{P}] = ^ E E P( X ) V ( X - * W) - M E ( 2 ) 

x=l x'=l x=l 

where the interaction potential is 

V(Ax) = - Ji (s(Ax + 1) + 5(Ax - 1)) + J 2 (5(Ax + 3) + 5(Ax - 3)) . (3) 

The Hamiltonian ([2]) can be rewritten in terms of the "unoccupancy" operator v(x) = 
1 — p(x) (0(x) = 1,0 for an empty and full site x respectively) 

H m\ = \ E E *(?M x - + 0* - ^o) E + L (y - y) ( 4 ) 

where V = J2 X V( x ) = 2(^2 — •A)- Note that Eq.([2]) in terms of p and Eq.(jl]) in terms of 
v have the same form for p = Vq/2. Moreover, the probability of the microstate {p(x)} for 
p = Vq/2 — Ap is the same as the probability of the "negative" of this microstate, {1 — p(x)}, 
for p = Vq/2 + Ap. Because of the above particle- hole symmetry the phase diagrams must 
be symmetric with respect to the symmetry axis p = J 2 — J\- 

We choose J\ as the energy unit and introduce dimensionless variables for any quantity 
X with dimension of energy as X* = Xj J 1; in particular 

T* = k B T/J u J* = J 2 /J 1: p* = p/J x . (5) 

B. The ground state 

The grand potential 

= - p L = -k B T In S = U -TS - pN (6) 

where p, U, S, N are pressure, internal energy, entropy and average number of particles re- 
spectively, reduces to the minimum of -ff[{p(x)}] for T = and fixed L. In the case of 
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periodic phases the bulk properties must be determined for L = In, where I is the period 
of density oscillations and n is integer. We consider u* = —p* = H*[{p(x)}]/ (In) for two 
homogeneous phases, one with all sites empty (gas) and the other one with all sites occupied 
(liquid) and for a periodic phase where three occupied neighboring sites are followed by / — 3 
empty sites with I > 6. For these phases we have 

empty lattice (gas) 

= < _2+3m1 periodic, I > 6 

J* — 1 — p* full occupancy (liquid) 

Two phases can coexist for thermodynamic states such that u* in these phases takes the 
same value. The (J*,p*) phase diagram for T* = is shown in Fig{TJ 

Note that for p* = —2/3 the u* of the periodic phase is independent of I if I > 6. This 
is because when in the empty lattice 3 neighboring cells become occupied, the associated 
change of H* is —2 — 3/i*. Because the interaction range is 3, for p* = —2/3 a triple of 
occupied cells can be separated from another triple of occupied cells by / — 3 empty cells for 
any I > 6. Such a state can be interpreted as a cluster fluid that can be stable, however, for a 
single value of the chemical potential, p* = —2/3. When 2 + 3p* > 0, the lowest value of u* 
corresponds to / = 6, and the periodic phase with period 6 is stable. The gas and periodic 
phases coexist for p* = —2/3. Because of this coexistence, arbitrary separation between 
the clusters, whose number is also arbitrary (but smaller than L/6), can be interpreted as 
arbitrarily small droplets (larger than 6 in the case of the periodic phase) of these phases. 
As a result, an arbitrary number of interfaces can be formed. This is possible when the 
surface tension between the gas and the periodic phases vanishes. 

Similarly, creation of a triple of empty sites in the fully occupied lattice leads to the 
change of H* which is AH* = —6 J* + 4 + 3/i*. At the coexistence between the fully 
occupied lattice and the periodic phase 6J* — 4 — 3p* = 0, hence the separation between 
the three empty neighboring cells (bubbles) can be arbitrary (but > 3). Again, such a 
state can be interpreted as a fluid of bubbles, or as a coexistence between the liquid and 
periodic phases in the case of vanishing surface tension. Note the similarity between this 
property of our model and the very low surface tension between water- or oil- rich phases 
and microemulsion. At T* = formation of the microemulsion is associated with vanishing 
surface tension in the lattice model for the water-oil-surfactant mixture 331 ]. 
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FIG. 1: Ground state of the considered model. The repulsion to attraction ratio J* and the 
chemical potential \x* are dimensionless (see ©). The coexistence lines are P* gas _ con( i = — 1 + J*, 
Moas-per = -2/3 and 

f^cond—per — 4/3 + 2J*. Schematic illustration of the three phases is shown 
in the insets inside the region of stability of each phase. 

III. MF APPROXIMATION 
A. Short background 

In the self-consistent MF approximation the Hamiltonian fl2]) is approximated by 



H MF [{p{x)}\ = J2 [ ~ (M*) + riftx) + -h(x)p( 

x 

where the mean-field acting on the site x is 

h(x) = -^V{x-x')p{x') 

x' 

and the MF average density satisfies the self-consistent equation 

e /9(/i(a;)+A0 



X 



(7) 



(8) 



p(x) 



The grand statistical sum 



n [«-* 



1 4. e P{h(x)+/j,) ■ 



§h(x)p(x) f 1 + e /3(h(x)+n) 



(9) 



(10) 



together with (jUJ) after some algebra leads to the grand potential of the form [34 1 

L L fl 

Xl =l x 2 = l L 

+ 5 Kr (xi - x 2 )fh(p(xi)) > -/i^p(x). 



(11) 



x=l 



where in the lattice models 



fh{p) = -k B Ts{p) = k B T pln(p) + (1 - p) ln(l - p) 



(12) 



Local minima of ( TTTj) satisfy Eq.(j9]) (see Ref.|34|). 

Eq.Q can be solved by iterations for different initial conditions. Stability regions of 
different phases and first-order transitions between them can be obtained by comparing 
u = Q/(ln) for different forms of p(x). In practice systems with the size I and PBC represent 
one period of the phases with the period I, and we have considered 6 < I < 50. 



B. Stability analysis 

Boundary of stability of the disordered phase can be found by analyzing the second 
derivative of (3Q with respect to the density profile p(x). The disordered fluid is stable as 
long as this derivative, 



C(X,X')= ^J" jr , (13) 



dp(x)dp(x') 

is positive definite. The disordered phase is at the boundary of stability when the smallest 
eigenvalue of ([TBI vanishes for p(x) = const. For interactions depending only on x — x' the 
quadratic part of (bilinear form) is diagonal in Fourier representation, and the eigenvalues 
of C are given by 



where 



In this model 



C(k) = FV{k) + -^—y ( 14 ) 



V(k) = ^V*(x)e ikx . (15) 



V(k) = -2 cos A; + 2 J* cos 3A;. (16) 

C(k) assumes the smallest value C{k^) for given T* and p* for k = corresponding to the 
minimum of V(k). We obtain 

/ if J* < 1/9, 

arccos \J otherwise, 

and 

' 2(J* - 1) if J* < 1/9, 

V(k b ) = i , \ 3/2 

-2J*fi±^-J otherwise. 
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FIG. 2: V(k) given by ([IB]). J* = 0.05, 1/9, 0.5, 1, 1.5 from the bottom to the top line on the left. 

Let us make a digression on the information on possible types of ordering that is included 
already in the form of V(k). For k = a macroscopic gas-liquid separation occurs (infinite 
wavelength 27r/0). For k > regions of the size 7r/k with excess number of particles, 
separated by depleted regions of the same size are formed. From the form of Q (see the 
first term in Eq. ffTT]) ) it follows that when in the disordered phase a density wave with 
the wavelength 2ir/k and the amplitude |0(A;)| is excited, then the energy increases by 
\<fi(k)\ 2 V {k) / Air . Note that for J* > 1 the energy of the homogeneously mixed phase is 
lower than the energy of the system separated into homogeneous gas and liquid phases, 
since V(0) = 2(J* — 1) > 0. Such separation is unfavourable energetically. For 1/9 < J* < 1 
we have V^O) < 0, but V(k) for k = assumes a maximum (Figj2j). It means that the 
gas-liquid separation is energetically favored compared to the homogeneously mixed sample, 
but the energy in the presence of the density wave with the wavelength 2irjk b is lower 
than in the case of macroscopic separation. Competition between the energy and entropy 
determines whether the periodic ordering or gas-liquid separation should occur in this case 
(1/9 < J* < 1) for given chemical potential and temperature. 

The boundary of stability of the disordered phase obtained from C{k b ) = 0, 

k B T* = -V(k b )p (l-p ), (17) 

represents the spinodal line of the gas-liquid separation when k b = (J* < 1/9), and the 
A-line marking the MF instability with respect to periodic ordering with the wavelength 
2n/kb for kj, > (J* > 1/9). It is also interesting to find the A-line and the spinodal in 
the (/i*,T*) phase space. From the form of the chemical potential for p = const, p* = 
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2( J* — l)p + T* ln[p/ (1 — p)], we obtain the boundary of stability of the homogeneous phase 

J* /l + 3J*\3/ 2 



p* = (J*-l)(l±g)+T*ln 



3J* 



where 



1 + 2: 



T* / 3J* 



3/2 



J* VI + 3J* 

The shapes of the spinodal and A-lines for various J* are shown in Fig]3] 



(18) 



(19) 




FIG. 3: Lines of instability (solid) of the homogeneous phase in the (p*,T*) variables for a range 



351 ] . The coexistence lines at T* = are shown as 



of J*. Similar behavior was obtained in Ref. 
dashed lines. 



In the case of attraction dominated system (J* < 1/9), the two branches of the spinodal 
line separate the region where Q assumes two minima for two different constant densities 
(low-T* side of the lines) from the region with one minimum on the high-T* side of the lines. 
At the left branch of the spinodal the minimum assumed for higher density disappears, 
whereas at the right branch the minimum assumed for the lower density disappears. The 
phase transition can occur on the low-T* side of the spinodal lines. This is usual behavior 
associated with the gas-liquid separation. 

For 1/9 < J* < 1 the branches of the line of instability intersect and form a loop for 
high T*. Inside the loop the two minima of Q corresponding to constant densities disappear 
(high-density minimum disappears to the left from the branch that corresponds to small \x for 
T* = 0, and the low-density minimum disappears to the right of the branch that corresponds 
to large \x for T* = 0). Thus, the homogeneous phase is unstable for any density inside the 
loop. 
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Note that when the loop is present, then for decreasing T* or increasing //* there may 
exist a sequence of phases: disordered - periodic - disordered. Such a sequence agrees 
with observations of the reentrant melting 0, 3^]. For the repulsion-dominated system 
(J* > 1) only the loop is present (FigfJ]), and coexistence between homogeneous phases is 
not expected. 

Our stability analysis is incomplete, because we did not study the boundary of stability 
of the periodic phase. 



C. Tricritical point 

We assume that when the continuous transition to the homogeneous phase with density 
p (p*,T*) is approached, the density in the periodic phase has the form p(x) = Po + Ap(x) 
with 

Ap(x) = 5p + (p cos(k b x) . (20) 

The continuous transition coincides with the A-line, and at the A-line 5p = = 0. The 
difference between the grand potential in the periodic and the homogeneous phases, AQ = 
AuL, is a function of 5p and <fi and Auj can be Taylor expanded for 5p — > and </> — > 0. When 
2ir/kb is not integer, in calculating Auj (see ffTTT) ) we choose n — > oo such that L m 2nn/kb 
and make the approximation 

cos(M m ~ — / cos z m dz. (21) 

The second derivative of Aoo with respect to <fi vanishes at the continuous transition, 
while the second derivative with respect to 5p is positive. From the extremum condition 
dAu/d5p = we obtain 

= JM 2 Q{ 4) (22) 

and 



where 



f3Au = a 2 <j) 2 + a 4 4 + 0(<f) (23) 



(3*V(k h )+A 2 (p ) 
a 2 (24J 
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a 4 




Mpo) 



(25) 



f3*V(0) + A 2 (p ) 

and A n (p) = d n (3fh(p)/dp n with f h given in Eq. ffT2l . The transition is continuous for a 4 > 0, 
and becomes first order at the tricritical point (TCP) given by a 2 = 0,4 = 0. We obtain for 
the density and temperature at the TCP the following expressions 



1 ± 



V(k b ) - V(0) 
V(k b ) + 3V(0) 



(26) 



-v{h)v(o) 



■ tcp 



(27) 



V(k b ) + 3V(0) 

Real positive solutions for p^ exist for 1/9 < J* < 1, ie. when the lines of instability 
intersect and form a loop. The T f * increases from for J* decreasing from 1, and the two 
TCPs merge at p = 1/2 and T* = 4/9 for J* = 1/9, i.e. when the disordered fluid starts 
to be unstable with respect to gas-liquid separation rather than with respect to periodic 
ordering. Thus, the two TCPs merged together change into the standard critical point for 
J* < 1/9. 



D. MF phase diagrams 

For very low temperatures the phase diagram is similar to the ground state (Fig{TJ). The 
gas and liquid phases are stable for very low and very high values of p* respectively. If 
J* < 1/3, the two fluid phases coexist. If J* > 1/3, the gas coexists with the periodic phase 
for some value of the chemical potential, and the periodic phase coexists with the liquid for 
some larger value of p*. The three phases coexist at a triple point (TP) for J* = 1/3 at 
T* = 0. 

Properties of the high T* part of the phase diagram can be deduced from the stability 
analysis. When J* < 1/9 only standard gas-liquid coexistence is expected. For J* > 1/9 
there is a continuous transition between the fluid and the periodic phase when T* > T^ cp and 
the transition becomes first order when T* < T^ cp . Since for high T* the periodic phase is 
stable for J* > 1/9 and for T* = it is absent for J* < 1/3, temperature at the TP increases 
from T* = to T* =4/9 for J* decreasing from 1/3 to 1/9. Note that for J* = 1/9 the TP 
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and the two TCPs merge, and this special point transforms to the standard critical point 
for J* < 1/9. 

For J* > 1/3 the TP and the coexistence of the homogeneous phases are absent. T 4 * p 
decreases for increasing J*, T^ cp = for J* = 1 and the two TCPs disappear for J* > 1. 
Thus, in the repulsion-dominated systems the high-T* transition to the periodic phase is 
only continuous. This seems to be inconsistent with the presence of the first-order transition 
between the periodic and the fluid phases at very low T*. This apparent inconsistency follows 
from the presence of two periodic phases for J* > 1. One of them is the same as the phase 
stable at T* = 0. It has large amplitude of density oscillations and the period I = 6. The 
other phase appears inside the loop of the A- line, has a period 2ir /&& and small amplitude of 
density oscillations. The two periodic phases coexist along the line which is a continuation of 
the coexistence line between the large-amplitude periodic phase and the homogeneous fluid, 
above the temperature at which this transition and the low-T* branch of the A-line (see 
Fig. 3) intersect. The coexistence between the two periodic phases terminates at a critical 
point, where the densities, amplitudes and periods of the two phases become the same. 

The precise locations of the first-order transition lines have been obtained by calculating 
the grand potential ffTTj) for the average densities that are self-consistent solutions of Eq.([9]). 
The method of determining the transition lines is shown schematically in FigJH We have 
chosen three characteristic values of J*, J* = 3, 1/3, 1/4. J* = 3 corresponds to the repulsion 
dominated system; no TP is present in this case. For J* = 1/3 the TP is present at T* = 0, 
and the repulsion-dominated system crosses over to the attraction-dominated system. For 
J* = 1/4 the TP is present at T* > 0, the two homogeneous phases coexist for low T*, and 
the periodic ordering occurs only at high temperatures. 

The (//*,T*) and (p,T*) phase diagrams for J* = 3 are shown in Figj5l The amplitudes 
of the two periodic phases along their coexistence line and for T* = 0.347 are shown in 
Figj6]as functions of T* and fi*. The density profiles in the two periodic phases for selected 
thermodynamic states are shown in FigJTl As far as we know, coexistence of two ordered 
phases with the same symmetry but diffrent degree of order has not been reported yet. 

The (fi*,T*) and (p,T*) phase diagrams J* = 1/3 are shown in Fig. [HI In this case 
there is a single periodic phase that coexists with the dilute and the dense homogeneous 
fluid for T* < T 4 * p and undergoes a continuous transition to the fluid for T* > T* cp . Similar 
phase behavior was obtained in Ref. 22j, where only one-dimensional density oscillations 
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FIG. 4: Illustration of the method used for obtaining the phase coexistence. In the bottom panel the 
left, central and right line corresponds to the gas, periodic and disordered liquid phases respectively. 



were assumed in the Landau-type and density-functional theories. 

The (/i*,T*) and (p,T*) phase diagrams for J* = 1/4 are shown in Fig. [9j In this case 
the two fluid phases coexist at T* < T t * ps 0.17063, while for T* > Tj* p each of them coexists 
with the periodic phase up to T t * p where the transition becomes continuous. For some range 
of temperature another periodic phase, with the same amplitude but with a larger period 
(hence lower density) becomes stable. This phase diagram is similar to the phase diagram 
obtained in Ref.[22] for very weak repulsion. Moreover, when electrostatic repulsion is added 
to the Landau functional, similar phase diagram is obtained 



IV. EXACT SOLUTIONS 



A. Transfer matrix and exact expressions 

Since the range of interactions is 3, we coarse-grain the lattice and introduce L/3 boxes 
consisting of three neighboring lattice sites, and labeled by integer 1 < k < L/3. The 
microstates in the k-th box are 

S(k) = (p(3k - 2), p(3k - 1), p(3k)). (28) 

There are 2 3 possible microstates in each box. We distinguish 4 states with the first site 
occupied and the remaining sites either occupied or empty, and denote such states by S±(k) = 

14 
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FIG. 5: MF phase diagram for J* = 3 in variables (p*,T*) (a) and (p, T*) (b). The symmetry 
axis is p* = 2 and p = 1/2 in (a) and (b) respectively. Only half of the phase diagram is shown 
because of the symmetry. Dash and solid lines represent continuous and first-order transitions. 
The dotted line is the A-line. The two-phase regions in (b) are shaded with different shades for 
different phase equilibria. The high-amplitude periodic phase coexists with gas (for p* < 2 or 
p < 1/2) or liquid (for p* > 2 or p > 1/2) for T* < T t * and with the low-amplitude periodic 
phase for T* > T t * . The line of continuous transitions from the homogeneous to the low-amplitude 
periodic phase terminates at the two first-order transition lines (one for p* < 2 or p < 1/2, the 
other one for p* > 2 or p > 1/2) at T^ p . The coexistence line between the two periodic phases 
terminates at the critical point with T* ~ 0.34713. 



(l,p(3/c — l),p(3&;)). Likewise, we denote states with the second site occupied by S2{k) = 
(p(3k — 2), 1, p(3k)), and with the third site occupied by Ss(k) = (p(3k — 2), p(3k — 1), 1). 
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FIG. 6: Amplitudes of the density profiles in the two periodic phases for J* = 3. (a) as a function 
of temperature along the coexistence line (the lines meet at T* ~ 0.34713) (b) as a function of \x* 
for T* = 0.347. 



The Hamiltonian of the system with PBC can be written in the form 

L/3 

H*[{p{x)}] = ^H* t [S{k),S{k + l)\ 



k=l 



where 



3fc 



X 



H;[S(k),S(k + l)]= [-p(x)p(x + l) + rp(x)p(x + 3)-f,*p( 

x=3k-2 

We introduce a 8 x 8 transfer matrix T with the matrix elements 

8 

T(S(k),S(k + l)) = e -P*Htm),s(k+i)] = J2 p ^(S(k))X i P l - 1 (S(k + 1)), 

8=1 



(29) 



(30) 



(31) 



where the eigenvalues of T are denoted by Aj such that |Aj| > |Aj + i|, the elements (S,i) 
of the matrix P transforming T to its eigenbasis are denoted by Pi{S), and the elements 
(i, S) of the matrix inverse to P by P i ~ 1 (5'). Note that T is not symmetric, hence pairs 
of complex-conjugate eigenvalues may occur. However, because T is a finite matrix with 
positive elements, from the Frobenius theorem it follows that the largest (in absolute value) 
eigenvalue is non-degenerate. 
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FIG. 7: Density profiles (a) in the coexisting high- and low-amplitude phases for T* = 0.3. The 
lines are shifted horizontally for clarity, (b) in the low-amplitude phase close to the continuous 
transition to the fluid at T* = 0.3 and (c) for T* = 1.3 and /i* = 1. The quasi-periodic structure 
with a period incommensurate with the lattice is obtained from a density profile with a large-period 
when 2-7r/fcb is noninteger. The lines connecting the results for integer x are to guide the eye. 



The grand statistical sum in terms of the transfer matrix takes the form 

5 = TrT L / 3 = >i /3 
t=i 



(32) 



and for the grand potential we obtain 



tt*/L = -p* = -T* 



i.nA 1 + Iln( 1+E (£ 

i=2 



X, \ L/3 



J -L->oc §" lnAl - 



(33) 



In the case of PBC the average density is independent of the position, p = (p(l)). From the 
definition of the average density and from Eqs. f l3~Tj) and fl32|) we obtain 



<P(1)> 



£&m £f=i ^ /3p i rl (S 1 (l))P i (S 1 (l)) 



8 x L/3 



£i=l \ 



il^oo ^(^(l))^^!^))- (34) 
Si(l) 
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FIG. 8: MF phase diagram for J* = 1/3 in variables (p*,T*) (a) and (p,T*) (b). The symmetry 
axis is fjL* = —2/3 and p = 1/2 in (a) and (b) respectively. Only half of the phase diagram is shown 
because of the symmetry. Dash and solid lines represent continuous and first-order transitions 
between the disordered fluid and the periodic phase. The dotted line is the A-line. The periodic 
phase in (a) is stable inside the loop (thick line). The two-phase regions in (b) are shaded. 

Let us consider the correlation function for the sites separated by a distance x = 3k + i, 
where k > and i — 0, 1, 2. Because the interaction range is 3, and the transfer matrix 
operates between triples of sites, we shall obtain different expression for (p(l)p(l + x)) for 
different i — 0, 1, 2. We introduce the notation 

G(3k + i) = (p(l)p(l + 3k + i)) - (p(l)) 2 . (35) 

From the definition of (p(x)p(x')) and from Eqs. fl3T|) and we obtain the asymptotic 
expression for L — > oo 



G(3k + ^) = Y: (v) A ^ B ^ 



where 2 = 0,1, 2, 



(36) 



(37) 
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FIG. 9: MF phase diagram for J* = 1/4 in variables (p*,T*) (a) and (p,T*) (b). The symmetry 
axis in (a) is //* = —3/4 and p = 1/2 in (b). Only half of the phase diagram is shown because of the 
symmetry. Dash and solid lines represent continuous and first-order transitions. The dotted line 
is the A-line. The coexisting phases in the two-phase regions in (b) are labeled by g for gas, d for 
dense fluid and p\, pi for the periodic phases with the smaller and the larger period respectively. 
The density range of stability of the large-period phase is within the thickness of the line. The 
periodic phases in (a) are stable inside the lense {jpi) and inside the loop (pi). 

B J U) = Y^Pn 1 (S j )P 1 (S j ), (38) 

and Sj is denned below Eq. (1281) . The asymptotic decay of correlations for k ^> 1 is deter- 
mined by the eigenvalue A 2 with the second largest absolute value. 
If A 2 is real, then for k 3> 1 and i = 0,1,2 we can write 

G(3k + i) = {sgn(X 2 )) k e- 3k ^A?B[% (39) 

where the correlation length is 

{ = 3/ in (A.). (40) 
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Note the qualitatively different behavior for A2 > and A2 < 0. For A2 < the correlation 
function changes sign when the separation between the particles increases by 3, in analogy 
with the density of the periodic phase in the ground state. The case A2 > corresponds 
to decay of correlations in the gas or liquid phases where no clusters consisting of three 
particles separated by 3 vacancies are formed. 

If A2 is complex, then A3 = A2, A^ = A± and = Bj 2 ^*. We introduce the notation 

A 2 = |A 2 |e iA , A {2) = |4 2 V Q1 , Bf = \Bf\e^, (41) 

and for k ^> 1 and i = 0,1,2 obtain the asymptotic expression 

G{3k + i) = A ie - 3k ^ cos (kX + 6$ (42) 

(2) (2) 

where Ai = 2\A\ ||-B^| and 9i = ot\ + 0i +i , i = 0, 1, 2. Similar expression was proposed in 
Ref. [3] for a 3d system. The structure factor obtained in experiments and theory [l, 7 



13- 



23j is also consistent with this form. In general, — n < A < tt, and 27r/A is noninteger. 
Except from A = ±7r (but in this case the imaginary part of A2 vanishes), the period of the 
exponentially damped oscillations is incommensurate with the lattice. This is similar to the 
results of the MF stability analysis and to the incommensurate density profiles obtained in 
MF for higher temperatures. 

When A = ±7r =p e with e <C 1, then we can write Eq. (T42|) in the equivalent form 

G(3k + i) = {-l) k e~ 3k/ tg{k, i) (43) 

with 

G(k, i) = Ai cos(6nk/w + (pi), (44) 
where the phase and the period of the amplitude modulations are <pi = —8iSgn(X) and 

w = H Tvil' ( 45 ) 

|A — s^n(Aj7r| 

The first factor in (143]) changes sign when the distance increases by 3. The last factor 
describes the modulated amplitude with the wavelength of modulations w ^> 6 if A — > ±7r. 

We have obtained A» and the matrix P numerically for different J*,n* and T* and the 
results are presented in the next section. 
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B. Results 



There are no phase transitions in a thermodynamic sense in one- dimensional systems. 
However, instead of a discontinuity, a rapid change of the density as a function of p* or 
p* can occur. Moreover, instead of long-range order and the associated periodic density, 
a short-range order with exponentially damped oscillatory decay of correlations with very 
large correlation length may exist. In order to verify if such pseudo phase transitions occur 
in this model, we calculate density and pressure for several values of J* for the range of p* 
and T* corresponding to the phase transitions obtained in MF. In the next step we examine 
the correlation functions. 

1. Thermodynamic properties (Equation of state) 

In Figs. irUlfTTI pf/i*) and p(p*) obtained from Eqs. f[3"3"j) and flM|) are shown for J* = 3 
and J* = 1/4 for 0.05 < T* < 1. By eliminating p* from Eqs. ([MD and ([MD we obtain the 
EOS, and present several isotherms p(p*) in FigsJT2] and [131 The chosen strengths of the 
repulsion to attraction ratio correspond to qualitatively different ground state and MF phase 
diagrams (see Figs. 1, [5] and [9]). Let us first discuss J* = 3. For low T* one can observe that 
although p*(p*) is a smooth function, its slope changes rapidly for the two values of p* that 
correspond to the phase transitions at T* = and in MF. In accordance with this behavior 
the density changes from nearly to 1 /2 in a very narrow range of p* and p* , remains nearly 
constant for large intervals of p* and p*, and again changes rapidly from p — 1/2 to p pa 1. 
Very large compressibility for p ^ 1/2 changes to very small compressibility for p pa 1/2. 
It is necessary to substantially increase the pressure in order to induce a slight increase of 
the density from p = 1/2, and a slight further increase of pressure is sufficient for a rapid 
compression to p pa 1. When T* increases from T* pa 0.1, the density changes from the gas 
density to 1/2 more and more gradually. For T* > 0.5 there are no abrupt changes of the 
slopes of the pip*) and p(p*) lines, but the curvature of these lines is significantly smaller 
than in the one-phase region of a simple fluid. This is because the repulsion between the 
particles at the distance 3 leads to a significant increase of pressure for random distribution 
of particles. On the other hand, small pressure for p < 1/2 when T* is low signals that in 
majority of states clusters made of at most 3 particles are separated by at least 3 empty 
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FIG. 10: p{fjL*) obtained from Eq.(|33]) (a) J* = 3; top line: T* = 0.1, bottom line: T* = 1 and (b) 
J* = 1/4; top line: T* = 0.05, bottom line: T* = 0.5. 




FIG. 11: pO*) obtained from Eq.flMD for J* = 3 and T* = 0.1,0.4,0.7, 1 (top to bottom line on 
the right) (a) and J* = 1/4 and T* = 0.005,0.05,0.1,0.15 (top to bottom line on the right) (b). 

sites. Similar behavior is observed for 1 < J* < 3 (Fig JTB"]) . but the range of //*, p* and T* 
for which p m 1/2 and remains nearly constant decreases with decreasing J*. For J* < 1 
the plateau at the p(p?) and p(p*) curves for p = 1/2 disappears. 

For J* = 1/4 we can see in Fig. [10] b) a rapid change of the slope of the p*(p*) line and 
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FIG. 12: EOS p(p*) isotherms obtained from Eqs.(|33D and ([Ml) for T* = 0.1,0.2,0.3,0.4,0.5 and 1 
(top to bottom line on the left) for J* = 3 (a) and J* = 1/4 (b). 
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FIG. 13: EOS p(p*) isotherms obtained from Eqs.flU and ([34]) for T* = 0.1. From the left to the 
right line J* = 0.1, 0.25, 1/3, 0.5, 0.75, 2.75, 3. 

in Fig. [11] the corresponding change of density from p ps to p ~ 1 when T* < 0.06. For 
T* > 0.15 the shapes of p*(p*),p(p*) and p(p*) (Fig jl~2]) resemble the corresponding curves 
in the single-phase simple fluid. We thus see a pseudo-transition between the gas and liquid 
phases for very low T*. We conclude that the thermodynamic properties show no signature 
of the weakly-ordered periodic phase previously found in MF (FigJHJ). By comparing FigsJl~2l 
a) and b) one can see the much lower pressure in this case than for the repulsion-dominated 
case of J* = 3. 
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2. Structure ( correlation function ) 

Our aim in this section is to discuss the exact results for the correlation function for 
J* = 3 and J* = 1/4, corresponding to qualitatively different ground state (Fig{TJ) and MF 
phase diagrams ( Figs \5\ and [H])- We particularly address the question for what parameters 
the periodic order occurs, and how the range and amplitude of the correlation function 
depends on /x*, T* and J*. 

For J* = 3 we obtain complex A2 for the considered region of (/i*,T*). In this case 
the correlation function is given in Eq. (|42p . and presented in FigsJT4l and [T5l In FigJT4l /1* 
corresponds to p ~ 1/2, where the periodic phase is predicted in MF, and in FigJTBI p* 
corresponds to p ~ (homogeneous gas in MF). 

a) 
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-0.2 
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90 96 102 108 114 120 126 
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FIG. 14: The correlation function G(x) for x = 3k + i with i = 0, 1, 2 (Eq.(@2])) for J* = 3, p* = 
and T* = 0.1 (inside the MF stability region of the periodic phase). Solid line and the circles 
(black), dashd line and the asterisks (red) and dotted line and the squares (blue) correspond to 
i = 0, 1,2 respectively. The bottom panel shows a small portion of the upper panel. 

The correlation length £ (Eq. (f40|) ) and the amplitude of the correlation function Aq (see 
below Eq. fj42l ) are shown in FigsJT6l and [TT1 respectively. For /i < —2/3 the wavenumber A 
is shown in Fig{18]a), and for /i > —2/3 the period w of modulations of the amplitude (see 
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FIG. 15: The correlation function G(x) for x = 3k + i with i = 0,1,2 (Eq.(g2])) for J* = 3, 
\i* = —0.7 and T* = 0.1 (outside the MF stability region of the periodic phase). Black (circle), red 
(asterisk) and blue (square) symbols correspond to % = 0,1,2 respectively. 

(gSD and OED) is shown in FigJT8]b). 
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FIG. 16: The correlation length £ (Eq.flU])) for J* = 3 as a function of T* (a) outside the MF stabil- 
ity region of the periodic phase. From the top to the bottom line /i* = —2/3, —0.7, —0.8, —0.9, —1 
and (b) inside the MF stability region of the periodic phase. From the bottom to the top line 
^* = -0.65, -0.6, -0.55, -0.5, -0.45 and 2. 

From Fig{16] a) it follows that for fi* < —2/3 the correlation length £ first increases 
slightly for decreasing T*, but starting form T* depending on /i* decreases rapidly to for 
T* decreasing to 0. Analogous behavior is predicted for /x* > 14/3 by the model symmetry. 
For —2/3 < /x* < 14/3 (stability of the periodic phase for T* = 0) the correlation length 
increases for decreasing T*. For given T* the correlation length increases with increasing 
jj,* when fi* < 2 and assumes a maximum for /j* = 2. The maximum of £ is very large for 
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T* < 0.15. For £ ~ 10 5 the range of the 'short-range order' is in fact macroscopic. For 
particles with a diameter lOnm the periodic arrangement persists to distances ~ 1mm. 
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FIG. 17: The amplitude .Ao of the correlation function (see Eq. (|42p and below) as a function of //*. 
Dash, solid and dotted lines correspond to T* = 0.1,0.2,0.3 respectively. J* = 3 (a) and J* = 1/4 
(b). 
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FIG. 18: (a) the wavenumber A of the correlation function (Eq. (|42|) ) . J* = 3 and /U* = 
-2/3,-0.7,-0.8,-0.9,-1 from the top to the bottom line (b) the period w (Eqs.(@3J)-(|lS])) of 
the amplitude modulations. J* = 3 and fi* = —0.65, —0.6, —0.55, —0.5, —0.45 from the bottom to 
the top line. 

The amplitude of the correlation function for T* < 0.1 increases sharply from a very small 
value for ji* < — 1 to ~ 0.2 for //* > —0.7 (see FigJTTh). The period w of the modulations 
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of the correlation function (Eqs.(l43|)- (145]) ) increases for decreasing T*, indicating more and 
more ordered structure (see Fig. [I8b). All these results confirm a qualitative change of 
the structure along the lines \i* = —2/3,14/3 for low T*. For —2/3 < ji* < 14/3 quasi 
long-range order with the very large correlation length and the amplitude that for low T* 
rapidly decreases at the boundaries of this region exists. From Figs {H] and [T7] we can see that 
the increasing correlation length and amplitude for increasing T* when — 1 < \i* < —2/3 
indicates a change from a less to a more ordered structure when T* increases. The rapid 
increase of the amplitude as a function of /i* for /i* —2/3 near T* corresponding to the 
maximum of £ resembles the transition between the periodic phases with weak and strong 
order found in MF (compare Figs. 17a andEb). When T* further increases, the properties 
of the correlation function change more gradually and the correlation length becomes short, 
in consistency with the continuous transition between the ordered and disordered phases 
found in MF for high T*. Despite the absence of the phase transition in the strict sense, we 
can see a change from a quasi-ordered periodic structure to the structure with much lower 
degree of order. 



a) 




T* 



FIG. 19: The wavenumber A (a) and the correlation length £ (Eq. (I40p ) (b) of the correlation 
function (Eq.(|42|)) as a function of T* for J* = 1/4. From the bottom to the top line in (a) and 
from the top to the bottom line in (b) p* = -0.75, -0.75 ± 0.01, -0.75 ± 0.02, . . . , -0.75 ± 0.07. 

For J* = 1/4 we obtain A2 that is a real number on the low-T* side of a line T* ross (/i*), 
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FIG. 20: The correlation function G(x) for x = 3k + i with i = 0, 1,2 (Eq.flM}) for J* = 1/4 and 
/i* = —3/4 for T* = 0.05 (a) and T* = 0.2 (b). Black (circle), red (asterisk) and blue (square) 
symbols correspond to i = 0, 1, 2 respectively. 

and a complex number for T* above this line. As a result, a monotonic decay of correlations 
is found for T* < T* ross (fi*), and an oscillatory decay sets in for T* > T* ross (fi*). The 
derivatives of the correlation length £ and the wavenumber A with respect to T* have a 
discontinuity when the imaginary part of A2 appears (see Fig JT9l) . The amplitude of the 
correlation function (FigJTTb) increases from a very small value for fi* < — 1 to ~ 0.3 for 
fi* = —0.75 for T* = 0.1. For larger T* the increase of the amplitude is more gradual. 
The correlation function shown in Fig. [20] confirms that for low T* the correlations decay 
monotonically, whereas for higher T* the oscillatory decay of correlations is present. The 
correlation length, however, is rather short, as shown in FigJT9l The monotonic decay of 
correlations for low T* , and oscillatory decay of correlations at higher T* together with a 
rather rapid increase of the amplitude of the correlation function from a very small to a 
rather large value for some range of /1* around //* = —3/4 bear some similarity to the MF 
phase diagram. However, due to the much shorter correlation length than for J* = 3, we 
can conclude that the weak periodic order (small amplitude of density oscillations) found in 
MF for relatively large T* does not resemble an ordered phase in Id. 





28 



V. MC SIMULATIONS IN CANONICAL ENSEMBLE 



In this section we present results of Monte Carlo simulations for the heat capacity. The 
simulations were carried out in the canonical ensemble. The basic step in the sampling 
is made as follows: Given the current configuration of the system one chooses at random 
with equal probability one of the occupied positions, x (with p(x) = 1), and one of the 
empty positions x' (with p(x') = 0), the trial configuration is then constructed by swapping 
the states between the positions x and x'. Considering the energies of the current and 



the trial configurations one applies the Metropolis criterion [39| to decide whether the trial 
configuration is accepted as the new configuration of the system or not. The heat capacity 
per particle, c v = (d(H/N)/dT) NjL , is computed using the fluctuation formula: 

* " [<" 2) " <* )2 1 ' (46) 

where the angular brackets indicate averages on the canonical ensemble. 

For J* = 3 the specific heat is shown in Figj2Tl The results for J* = 0.25 are shown 
in Fig. |22j For both values of J* we can observe the presence of a peak at low density. 
On cooling the system the height of this peak increases, and the density where it appears 
is shifted to lower values. This maximum in the heat capacity at low temperature can be 
explained as an effect of the equilibrium between isolated particles and clusters of several 
particles. Given the Hamiltonian of the model, these clusters are likely to be triples in the 
case of J* = 3. At low temperature and low density the loss of entropy due to the formation 
of clusters is compensated by the energetic effect due to the attractive interaction between 
the nearest neighbors. 

There are, however, significant differences between the heat capacities curves for J* = 0.25 
and J* = 3.0, especially for p ~ 1/2. For J* = 3.0, at low temperature, we can observe 
a basin around p = 1/2, and a narrow peak centered also at p = 1/2. Focusing in the 
region p < 1/2, the ground state configurations are formed by triples of occupied positions. 
Each triple of ocupied cells is separated at least by three empty positions from another 
triple. Since the triple-triple interaction is repulsive at short distances, the system does not 
show any trend to exhibit a pseudo phase separation to form large regions of occupied and 
empty positions, and therefore the small energy fluctuations lead to small values of the heat 
capacity. Notice, however that for J* < 1/3 and low temperature the dominant attractive 
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interactions lead to a condensation of particles in large clusters of occupied cells. In this 
case neither the basin at low T* nor the peak near p = 1/2 is present. 

The peak of c v (p) at p ~ 1/2 for J* = 3 can be interpreted as a signature of a pseudo- 
phase transition between an ordered (or quasi-ordered) phase (periodic phase with I = 6 
for p = 1/2) and a high temperature disordered phase. Notice that as one approaches 
p = 1/2 the degeneracy of the ground state reduces sharply, then we can describe this peak 
as produced by the competition between the ground state (with very low entropy values 
when p — > 1/2) and disordered states (with higher values of energy and entropy). This 
ordered pseudo-phase lies between the fluid of small droplets (p < 1/2, T — > 0) and the 
fluid of bubbles (p > 1/2, T — > 0). In spite of the lack of real phase transitions for one- 
dimensional models with short range interactions, the periodic pseudo phase with density 
p = 1/2 resembles to some extent the low-density crystalline phases that appear in core- 



softened models 



4Q| 



We conclude that measurements of c v in systems with competing interactions can give 
information on the formation and properties of clusters for very small densities, and on 
formation of phases (or pseudo-phases) with periodically ordered clusters for higher densities. 
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FIG. 21: The specific heat per particle (in ks units) as a function of density (dimensionless) for 
J* = 3 T* = 0.25, 0.5, 0.75, 1 with L = 1200 . 
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FIG. 22: The specific heat (in fee units) as a function of density (dimensionless) for J* = 1/4 and 
T* = 0.1, 0.25, 0.5, 0.75, 1 with L = 840. 

VI. SUMMARY AND DISCUSSION 



We have developed a generic model for self-assembly in systems with competing interac- 
tions. The Id version of the model was solved in MF approximation and exactly in the grand 
canonical ensemble for the whole range of the repulsion to attraction ratio J*. In addition, 
MC simulations have been performed in the canonical ensemble. Previously lattice mode 



with competing interactions were considered in the context of magnetic systems |31l. I41L 142] . 
The studies, however, did not focus on the role of the external magnetic field that in the 
fluid version of the lattice models corresponds to the chemical potential fi*. In the context of 
fluids the chemical potential plays a crucial role and allows to obtain the EOS and structure 
for dense and diluted fluids. We have found interesting and counterintuitive results even in 
Id case. 

The ground state (T* = 0) shows a sequence of phases fluid (gas) - periodic - fluid (liquid) 



for increasing /j,* when J* > 1/3. T 



experimentaly in several systems 



lis behavior agrees with the reentrant melting observed 
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371 ] . For weaker repulsions only the two fluid phases are 



present for T* = 0. A peculiar property of the MF solutions is the existence of the periodic 
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phase for a range of //* that is broader for intermediate temperatures than at T* = 0. In the 
case of 1/9 < J* < 1/3 th periodic phase appears for some range of T*, even though it is 
absent for T* = 0. Usually, the increase of T* leads to less ordered structures, and the MF 
result is counterintuitive. In particular, the MF phase diagrams show that for decreasing 
temperatures at constant /i* a sequence of phases fluid-periodic-fluid appears, i.e. we find 
reentrant melting. The high- T* periodic phase is characterized by a small amplitude of 
density oscillations and a period that is incommensurate with the lattice. 

There are no phase transitions in Id systems. In order to find if the MF phase behavior 
is associated with a qualitative change of mechanical and structural properties, we have 
analyzed the exact results for the EOS and correlation function. We have found that for 
J* > 1/3 the shapes of the EOS isotherms and the very large correlation length indicate 
pseudo phase transitions to the periodic phase. On the other hand, the high-T* weakly 
ordered phase found in MF is only reflected in the qualitative change of properties of the 
correlation function. For ji* that in the ground state corresponds to the gas or liquid, 
we observe that the correlation length assumes a maximum for some finite T*, indicating 
increasing order for increasing T*. Moreover, for such T* the amplitude of the correlation 
function changes from a very small value to a much larger value for a narrow range of fi*, 
and stays large for the range of ji* similar to the stability region of the periodic phase found 
in MF. The most amazing behavior shows the correlation function for 1/9 < J* < 1/3. 
It crosses over from a monotonic decay at low T* to an oscillatory decay for higher T* 
that is, however, much smaller than the temperature correponding to the appearance of 
the periodic phase in MF. The unusual appearence of the periodic short-range order at 
intermediate temperatures that for weak repulsion is neither energetically nor entropically 
favored is the most surprising and interesting exact result of this work. It means that 
for competing interactions the inhomogeneities may appear as a compromise between the 
macroscopic phase separation that is energetically favorable but entropically unfavorable, 
and the disordered structure that is favorable entropically and unfavorable energetically. 
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